Exclusive pp — > pp7r + ir reaction: 
from the threshold to LHC 

P. Lebiedowicz^B and A. Szczurek 2 ' 1 '0 

1 Institute of Nuclear Physics PAN, PL-31-342 Cracow, Poland 
2 University of Rzeszow, PL-35-959 Rzeszow, Poland 

Abstract 

We evaluate differential distributions for the four-body pp — ► ppn + n~ reaction which constitutes 
a irreducible background to three-body processes pp — ► ppM, where M are a broad resonances in 
the 7T+7T- channel, e.g. M = a, p°, / (980), / 2 (1275), / (1500). We include both double-diffractive 
contribution (both pomeron and reggeon exchanges) as well as the pion-pion rescattering con- 
tribution. The first process dominates at higher energies and small pion-pion invariant masses 
while the second becomes important at lower energies and higher pion-pion invariant masses. The 
amplitude(s) is(are) calculated in the Regge approach. We compare our results with measured 
cross sections for the ISR experiments at CERN. We make predictions for future experiments at 
PANDA, RHIC, Tevatron and LHC energies. Differential distributions in effective two-pion mass, 
pion rapidities and transverse momenta of pions are presented. The two-dimensional distribution 
in (y w + , y n - ) is particularly interesting. The higher the incident energy, the higher preference for 
the same-hemisphere emission of pions. The processes considered constitute a sizeable contribution 
to the total nucleon-nucleon cross section as well as to pion inclusive cross section. 
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I. INTRODUCTION 



Diffractive processes although very difficult from the point of view of perturbative QCD 
are very attractive from the general point of view of the reaction mechanism. There are 
several classes of diffractive-type processes [l[ in high-energy nucleon-nucleon collisions such 
as: 

• (a) elastic scattering, 

• (b) single- diffractive excitation of one of the nucleons, 

• (c) double-diffractive excitation of both participating nucleons, 

• (d) central (double)-diffractive production of a simple final state. 

The energy dependence of the first three types of the reaction was measured and can 
be nicely described in a somewhat academic two-state (but fullfilling unitarity) Good- 
Walker model The last case was not studied in too much detail either experimentally or 
theoretically. At not too high energies the dominant diffractive final state is the 7r + 7r~ and 
7r°7r° continuum. The multi-pion and KK continuum is expected to be smaller. 

There is recently a growing interest in understanding exclusive three-body reactions 
pp — > ppM at high energies, where the meson (resonance) M is produced in the central 
rapidity region. Many of the resonances decay into tttt and/or KK channels. The represen- 
tative examples are: M = a, p°, / (980), 4>, /2(1275), /o(1500), Xc(0 + ). It is clear that these 
resonances are seen (or will be seen) "on" the background of a tttt or KK continuum 1 . 
Therefore a good understanding of the continuum seems indispensable. In the present anal- 
ysis we concentrate on the 7r + 7r~ channel. Similar analysis can be done for 7r°7r° exclusive 
production. 

At larger energies two-pomeron exchange mechanism dominates in central production 
(see |l| and references therein). In calculating the amplitude related to double diffractive 
mechanism for pp — > pp7T + 7T~ we follow the general rules of Pumplin and Henyey [ij (for 
early rough estimates see also Ref. [E|). 

At lower energies subleading reggeons must be included in addition to the pomeron 
exchanges. We include a new mechanism relevant at lower energies (FAIR, J-PARC) relying 
on the exchange of two pion. We shall call this mechanism pion-pion rescattering for brevity. 

We discuss interplay of all the mechanisms in a quite rich four-body phase space. 



II. THE ttN ELASTIC CROSS SECTION 

At low energies, the total cross sections for ir + p and ir~p show a significant energy- 
dependent asymmetry defined as: 

The total cross section tests, via optical theorem, only imaginary part of the scattering 
amplitude. In our case of the 2 — > 4 reactions 2 we should use rather full scattering amplitude. 



1 In general, the resonance and continuum contributions may interfere. This may produce even a dip. A 
good example is / (980) production (see @ 5 0])- 

2 2 — > 4 reaction denotes a type of the reaction with two initial and four final particles. 
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In contrast to the total cross section the elastic scattering cross sections for ir + p and tt p 
show at low energies rather small asymmetry defined as: 



\a^ p (W) - a n , P (W)\ 
A*f(W) = 1 el + { ' el { n . (2.2) 
a^ p (W) + a^ p (W) 

A reliably model should explain such details of the interaction. 

Therefore to fix parameters of our double- diffract ive model we consider first elastic pion- 
proton scattering. The amplitude for the elastic scattering of pions on nucleons is written 
in the simplified Regge-like form: 

/ s \ ajp(t)-l / b jp 
M^± p ^± p {s,t) = isCpl— J expf-^t 

/ \ am(t)-l / tdM 

+ («/ + i)sC, exp [-f-t 

± (a p -i)sC p l^-j exp i^tj , (2.3) 

where a/ = -0.860895 and a p = -1.16158. The strength parameters Cjp, Cf, C p are taken 
from the Donnachie-Landshoff model |9l for total cross section: 



Cjp = 13.63 mb, C f = 31.79 mb, C p = 4.23 mb . (2.4) 

This means that our effective phenomenological model describes the available total cross 
sections. The pomeron and reggeon trajectories determined from elastic and total cross 
sections are given in the form (a,(t) = aij(0) + a^i): 

ajp(t) = 1.088 + 0.25t, a m (t) = 0.5475 + 0.93t . (2.5) 

The values of the intercept ajp(0) and ajR(O) are also taken from the Donnachie-Landshoff 
model [9] for consistency. The effective slope parameter can be written as 

B eff = B nN (W wN )=B + 2a[\n(j-) . (2.6) 

We take a[ = 0.25/0.93 for pomeron and reggeon exchanges, respectively. The slope param- 
eter -Bttat; taken the same for the pomeron and reggeons, must be fitted to the data. From 
the fit to the data [s| we find B = 5.5 GeV -2 . The effective slope observed in ^-distributions 
is of course much larger (B eff = 7-10 GeV" 2 for P lab = 3-200 GeV @). 

The differential elastic cross section is expressed with the help of the scattering amplitude 

as: 

:\M(s,t)f. (2.7) 



dt lQns 2 

The differential distributions da e i/dt for both ir + p and 7t~p elastic scattering for three 
incident-beam momenta of = 5 GeV, = 50 GeV and P^ = 200 GeV are shown in 
FigfTJ Under a detailed inspection one can observe that the local slope parameter 

d . ( do, 
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FIG. 1: Differential distributions for ir + p (left) and n~p (right) elastic scattering for different 
energies calculated with the amplitude (|3.2p and parameters as given in the text. In this calculation 
the slope parameter was taken as B w n = 5.5 GeV -2 (dashed line). A fit to the world ttN elastic 
scattering data suggest that the pomeron and reggeon slopes may be slightly different. The solid 
line shows such a result. The details are explained when discussing Figj2j The experimental data 
are taken from Ref.Q]. 

is t-dependent and is slightly larger for n~p than for 7i + p. Such an effect was observed exper- 
imentally in Ref.Q]. The local slope decreases with increasing t. A rather good description 
of experimental da e i/dt is achieved. 

Our one-parameter (B^n) model here is consistent with the simple Donnachie-Landshoff 
model for total cross section A more refined model should include absorption effects 
due to pion-nucleon rescatterings. The analysis of absorption effects clearly goes beyond the 
scope of the present paper. Our model sufficiently well describes the ttN data and includes 
absorption effects in an effective way. This has advantage for the pp — > ppim reaction 
discussed in the present paper where the irN absorption effects do not need to be included 
explicitly. This considerably simplifies the calculation for the 2— >4 reaction and actually 
this makes the calculation of the 2— »4 reaction feasible. 

Before we shall go to the pp —>■ ppir + 7r~ reaction, we have to discuss the parameters of the 
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ttN interaction. The strenght parameters of the pomeron and reggeon couplings are taken 
from the Donnachie-Landshoff analysis of the total cross section in several hadronic reactions 
[s| as discussed above. The only free parameters - the slope parameters, are adjusted to the 
elastic 7t + p and 7r~p scattering. With Bp = 5.5 GeV" 2 and B R = 4 GeV -2 we nicely describe 
the existing experimental data for np scattering as can be seen from Fig [2] (solid lines). The 
long-dashed lines show pomeron (IP) and reggeon (M) contributions and the short-dashed 
lines their interference term. In the Regge approach, high energy cross section is dominated 
by pomeron exchange. The reggeon exchange dominates in the resonance region. There is 
a region of energies where the interference term dominates. This is very different than for 
the total cross section where the cross section is just a sum of the pomeron and reggeon 
terms. We get a nice description of the data for y/s > 2.5 GeV. The region below contains 
resonances and is therefore very difficult for modeling. 
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FIG. 2: The integrated cross section for ttN elastic scattering. The experimental data are taken 
from Ref. 171 ]. 



Having fixed the parameters we can proceed to our four-body pp — > ppir + ir reaction. 



III. CENTRAL DOUBLE DIFFRACTIVE CONTRIBUTION 



Pa 




M 

7V + (p 3 ) 





FIG. 3: A sketch of the dominant mechanisms of exclusive production of tt + tt pairs at high 



energies. 



The general situation is sketched in FigfSl The corresponding amplitude for the pp 



pp7r + 7r process (with four-momenta p a + pb — > p\ + p 2 + £>3 + P4) can be written as 

^fl^Hwnr = Ml3 (t 1 ,S 13 )F(t a ) j - F(t a )M 2i (t 2 , S 24 ) 

+ M 14 (t 1 ,5 14 )F(t & ) 1 2 F(t 6 )M 23 (t 2 ,g 23 ) , (3.1) 

where Mjfc denotes "interaction" between nucleon z = 1 (forward nucleon) or i = 2 (backward 
nucleon) and one of the two pions k = 3 (7r + ), = 4 (7r~). In the Regge phenomenology 
they can be written as 

\ Ojp(ti)-l / R 



M ik (U, s ik ) = i Cp ( — J exp ( — t, 

+ (a/ + i) s ifc C/ ^^^j exp 



± (a p - i) s ifc C p I — J exp ( — *i 1 . (3.2) 

Above Sik — Wf k , where Wik is the center-of-mass energy in the (i, k) subsystems. The 
third term is with the sign plus if k = 3 and with the sign minus if k = 4. The normalization 
constants {Cp, Cf, C p ) can be estimated from the fit to the total 71N cross section (12.41) . The 
values of the Regge trajectories (I2.5P are also taken from the Donnachie-Landshoff model 

0. The first term describes exchange of the leading (pomeron) trajectory while the next 
terms describe the subleading reggeon exchanges. At high irN subsystem energies Wik > 20 
GeV only the pomeron exchange survive. 

The extra form factors F(t a ) and F(tj,) "correct" for off-shellness of the intermediate 
pions in the middle of the diagrams shown in FigOJ In the following they are parametrized 
as 

F( tl , 2 ) = exp (^f^) , (3.3) 

1. e. normalized to unity on the pion- mass-shell. In the following for brevity we shall use 
notation t%,2 which means t\ or t 2 . In general, the parameter h. ff,E is not known but in 
principle could be fitted to the (normalized) experimental data. From our general experience 
in hadronic physics we expect h. ff,E ~ 1 GeV. How to extract h. ff,E will be discussed in 
the result section. 

The parametrization j^] can be used only for Wik > 2 — 3 GeV. Bellow Wik = 2 GeV 
resonances in ttN subsystems are present. In principle, their contribution could and should 
be included explicitly 3 4 . The amplitude (13.11) with f)3.2p is used to calculate the correspond- 
ing cross section including limitations of the four-body phase-space. To exclude resonance 
regions we shall "correct" the parametrization ( 13. lj) with ( 13.21) by multiplying by a purely 
phenomenological smooth cut-off correction factor: 

wN exp ( w - w ° ) 

fcontiWik) = - | — ^ w-WP j ■ ( 3 ' 4 ) 



3 The higher the center-of-mass energy the smaller the relative resonance contribution. 

4 In the standard terminology the resonances belong to single-diffractive contribution to be distinguished 
from double-diffractive contribution discussed here. 
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The parameter Wo gives the position of the cut and parameter a describes how sharp is the 
cut off. The first parameter can have a significant influence on the results. We shall take 
Wo = 2 — 3 GeV and a = 0.1 - 0.5 GeV. For large energies f££ t {Wik) ~ 1 and close to 
kinematical threshold W ik = m n + M N : f££ t (Wik) ~ 0. 

IV. PION-PION RESCATTERING 

For the pp — > pp-K + 7i^ or pp — > pp-n^n^ reactions there is another type of semi-diffractive 
contribution shown in FigJH 




FIG. 4: A sketch of the high-energy pion-pion rescattering mechanisms. 



Similarly as for the pp — > A^A^/ (1500) reaction (see [131 ]) the amplitude squared - averaged 
over the initial and summed over the final state - for these processes can be written as: 



\M\< 



(E a + m) {E 1 + m 



Pa 



Pi 



2p a ■ Pi 



a 2 



(*1 



(£ a + m) 2 (Ex + m) 2 ( J E a + m)(E 1 + m 

F^NN^l) X |-^7r0*7r°*^7r+7r-(S34,^0;^l,^2)| 2 X 



fl 2 



(E b + m) (£ 2 + m) 



P^ 



+ 



Pi 



(*2 - n> 
2pb • p 2 



p2 
2 \ 2 r rcNN 



X 2 



(£ 6 + m) 2 (E 2 + m) 2 (£ 6 + m)(£ 2 + m) 



x 2 . 
(4.1) 



In the formula above m is the mass of the nucleon, E a , E^ and Ei, E 2 are energies of initial 
and outgoing nucleons, p Q ,P& and Pi,P2 are corresponding three-momenta and m % is the 
pion mass. The factor g n jy n is the familiar pion nucleon coupling constant and is relatively 
well known Q = 13.5 - 14.6). In our calculations the coupling constants are taken 

as gl NN /4ir = 13.5. 



0*^0* 



At high-energies the pion-pion scattering amplitude of the subprocess n n 
with virtual initial pions can be written, similarly as for ttN scattering: 



7T 7T 



(^34, to] tl, t 2 ) 



1J s 34 



so 



"«(*o)-l 



exp 



t J ^0.(^)^0.^2) . 

(4.2) 
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We have parametrized the to dependence in the exponential form. The slope parameter 
is not well know, however it may be expected to be B„ n ~ 4-6 GeV -2 . In the formula 
above i 7 ^* (^1,2) are extra correction factors due to off-shellness of initial pions. We use 
exponential form factors of the type (13.31) . To exclude resonance regions we "correct" the 
Regge parametrization (14. 2 p by multiplying by a factor f^ t (Wu) (as in (13.41) ). 

In the case of central production of pion pairs not far from the threshold rather large 
transverse momenta squared ti and t 2 are involved and one has to include non-point-like and 
off-shellness nature of the particles involved in corresponding vertices. This is incorporated 
via -Fttaw^i) or F^^^ih) vertex form factors. In the meson exchange approach [15( they 
are parameterized in the monopole form as 

A 2 _ irr,^ 

F*NN{tw) = ■ (4.3) 

Typical values of the form factor parameters are A = 1.2-1.4 GeV [151 ]. however the Gottfried 



Sum Rule violation prefers smaller A 0.8 GeV [161 . 



V. THE DIFFERENTIAL CROSS SECTION 

The differential cross section for the 2^4 reaction is given as 



d 3 pi d 3 p2 d 3 p 3 d 3 p4 



do = — \M.\ 2 (2ir) 4: 5 4: (p a + ph — pi — p 2 — Pa — Pa)- — — - — -7T- — -, — r4— — ; — (5.1) 
2s 1 lv ' yF F F F Fi y ' (2n) 3 2E 1 (2tt) 3 2£ 2 (2n) 3 2E 3 (27r) 3 2£ 4 

This can be written in a useful form: 

1 



da = —\M\ 2 5 A (p a +p b -p 1 -p 2 -p3-p4)- 



2« ^" ^ ^ ^ ™ ^ y (27r) 8 2 4 
x (dy 1 p lt dp lt d(j) 1 )(dy2P2tdp2td4>2)(dy 3 d 2 p 3t )(dy4d 2 p4t) ■ (5.2) 

This can be further simplified: 

da = ^-\Mf*5{E a + E b -E 1 -E 2 -E 3 - E A )5 3 (p lz + p 2z + p 3z + 1 ' 



x {dyipitdp lt d(j)i)(dy2P2tdp2td(t)2)dy 3 dy A d 2 p m . (5.3) 
Above we have introduced an auxiliary quantity: 

Pm = V-it ~ P4i • (5-4) 

We choose transverse momenta of the outgoing nucleons (pu,P2t,P3t,P4t), azimuthal angles 
of outgoing nucleons (0i, 2 ) and rapidity of the pions (y 3 , y 4 ) as independent kinematically 
complete variables. Then the cross section can be calculated as: 



, rr-\( 1 ; / \ | \M {Pit , 01 , P2t , 02, Vz , 2/4, Pra , 0m) 1 1 

d(J = >J (Plt,(Pl,P2t,<P2,y3,y4,Pm,(Pm)\k / . A = 777^7^7 

2y/s(s -Am 2 ) [Ztt)* 2* 

x pudp lt d(j)iP2tdp2td(t)2^dy 3 dyAd 2 p m , (5.5) 



8 



where the 5 functions have been totally eliminated and k denotes symbolically discrete 
solutions of the set of equations for energy and momentum conservation: 



~P3z - 



Piz 



™lt + Viz + V m lt + Plz 



Plz + P2z 



(5.6) 



where m\ t and m^t are transverse masses of outgoing nucleons. The solutions of Eq. (l5.6j) 
depend on the values of integration variables: p\ z = Pi z (pit,P2t,P3t,P4t, <fii, <f>2, V3, Va) and 

P2z = P2ziPlt,P2t,P3t,P4tAlA2,y3,yi)- 

In Eq. ( 15. 51) an extra Jacobian of the transformation (2/1, 2/2) ~~ > (piz,P2z) has appeared: 



Jk 



Piz(k) 



P2z(k) 



V m lt+Piz{k) 2 ^m 2 2t +p 2z (ky 



(5.7) 



In the limit of high energies and central production, i.e. p\ z ^> (very forward nucleonl), 
—P2z (very backward nucleon2) the Jacobian becomes a constant J — > |. 

To calculate the total cross section one has to calculate the 8-dimensional integral (see 
Eq. (15.51) ) numerically. This requires some care. 

In the next section we shall show our predictions for several differential distributions in 
different variables. 



VI. RESULTS 



Before we go to our four-body reaction let us focus for a moment on 7r°7r° — > 7i + n on-shell 



7T 



scattering. In Figl5]we show the total (angle-integrated) cross section for the 7rV -> 7T 
process. We include both the pion-pion rescattering contribution obtained from partial wave 
analysis 0] as well as contribution from the Regge phenomenology at higher energies. The 
parameters of the Regge amplitude for the irn — > tttt scattering were obtained in Ref.[Io| 
from different isospin combinations of nucleon-(anti)nucleon, and pion-nucleon scattering 
assuming Regge factorization. For our case of 7r°7r° — > 7r + 7r~ reaction only the p-reggeon 
exchange is relevant. We show predictions for the Regge contribution for corrected (Wo = 
1.5, 2 GeV and a = 0.2 GeV in Eq. (l3.4p ) extrapolations to low energies and for different 
values of the slope parameter B W7r = 4 GeV -2 (dotted lines), B nn = 5 GeV -2 (dashed line) 
and B n7T = 6 GeV -2 (solid lines). A relatively good matching is achieved without extra 
fitting the model parameters. In the following we shall focus on the higher-M^ Regge 
component which dominates at higher energies (see a next figure). 

In Figj6] we present the total cross section for the pp —> ppTT + n~ reaction, i.e. the cross 
section integrated over full phase space, as a function of the center-of-mass energy. We show 
theoretical predictions from the models calculations with A = 0.8 GeV and = 1 GeV 2 

(solid lines) and &off,E = 0-5 GeV 2 (dashed lines). The bottom dotted line was obtained 
with A = 0.8 GeV and A 2 off E = 0.5 GeV 2 while the top dotted line with A = 1.4 GeV and 
A-offE = ^ GeV 2 . Details of the low-M^ rescattering contribution can be found in Ref.0]. 
The search for a double pomeron exchange mechanism contribution leads to an upper limits 



of ~ 20/ib (for M n7T sC 0.7 GeV) [18j, (49 ± 5.5)/xb [19|, (30 ± ll]ub |20| and (44 ± 15)/xb 
2l| . The experimental value of the cross section taken from [20] was obtained for M p n 
> 2 GeV and no limitation on M W7V ; reduces however to 9 pb for M nn ^0.6 GeV [20] . 
For comparison we show the full cross sections for the pp — ► pp7r + 7r~ reaction (filled black 
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FIG. 5: The angle-integrated cross section for the reaction 7r°7r° — > 7r + 7r _ . We present contributions 
obtained from partial wave analysis [f| and Regge phenomenology 10] for corrected (Wq = 1.5, 2 
GeV and a = 0.2 GeV in Eq. (|3.4|) ) extrapolations to low energies. 



circles) from Ref. 22| and for the pp — > ppTT + 7r~ reaction (filled blue triangles) from Ref. 23 
which are more than 1 mb for (2.5 < y/s < 10) GeV 5 . Clearly for low energies (y/s < 20 
GeV) neither exclusive double diffraction nor pion-pion rescattering constitute the dominant 
mechanism. Here the production of single and double resonances is the dominant mechanism 
(see e.g. (sj). The mechanism of the resonant production is rather complicated and will not 
be discussed in the present analysis. 




io z 

\fs (GeV) 



FIG. 6: Cross section for the pp — > ppTr + ir~ reaction integrated over phase space as a function of the 
center-of-mass energy. We compare the pion-pion rescattering and double-diffractive contributions 
with the experimental data (open symbols represent DPE contribution from Refs. 18|, 1 19 . 
and filled symbols show the cross sections for the pp — > ppn + TT~ reaction (black circles) from Ref. 22| 
and the pp — ► ppir + ir~ reaction (blue triangles) from Ref. [3] ). The theoretical uncertainties for 
these contributions are shown in addition. 



The results depend on the value of the nonperturbative, a priori unknown parameter of 
the form factor responsible for off-shell effects. In Table [J we have collected integrated cross 



This is a significant contribution to the total pp cross section. 
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TABLE I: Full-phase-space integrated cross section (in /xb) for exclusive double diffractive tt + tt~ 
production at selected center-of-mass energies and different values of the off-shell-form factor pa- 
rameters. Here Wo = 2 GeV and a = 0.2 GeV in Eq. (|3.4p . No absorption effects were included 
explicitly. 



^1,2) 


Kff (GeV 2 ) 


W = 5.5 GeV 


W = 200 GeV 


W = 1960 GeV 


W = 14 TeV 


ew{(h,2-ml)/A 2 off ^ 


0.5 


0.1 


50.3 


96.4 


179.1 




1 


0.6 


146.2 


287.2 


535.2 


( A o//,M-^)/( A o//,M-*l,2) 


0.5 


0.02 


18.9 


35.6 


66 




1 


0.18 


64.6 


125.2 


232.8 


(( A o//,D-0/( A o//,D-M)' 


0.5 


0.31 


83.6 


164.2 


306.5 




1 


1.15 


217.5 


437.9 


822.4 



sections for selected energies and different values of model parameters. We show how the 
uncertainties of the form factor parameters affect our final results. 

In FigJ7|we show predictions for different values of the parameter A 2 jj E = 0.5 GeV 2 (lower 
lines), A o//,_e; = 1 GeV 2 (upper lines) and for naive (dashed lines) and corrected (solid lines 
with Wq = 2 GeV and a = 0.2 GeV) extrapolations to low energies. The experimental 
cuts on the rapidity of the pions are included when comparing our results with existing 
experimental data. Although not all the data are in good agreement with the predictions, 
their general trend follows the theoretical expectations. No absorption effects were included 
in this calculation. In general, the higher energy the higher absorption effects. The bare 
cross section rises with energy. The absorption corrections are expected to lower or even 
stop the rise. Consistent inclusion of absorption effects is rather difficult and will not be 
studied here. 
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FIG. 7: Cross section for the otj 
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-» jOj07r + 7r reaction integrated over phase space with cuts relevant 
25j, |26J]. The experimental value from [26J was obtained for the 



> 2. We show results for different values of the parameter A 2 



ff,E 



0.5 



for a given experiments 
different cut Ay = \y p — y 

GeV 2 (lower lines), A 2 ^ E = 1 GeV 2 (upper lines) and for the naive (dashed lines) and corrected 
(solid lines with Wq = 2 GeV and a = 0.2 GeV) extrapolations to low energies 



The distribution in the (y 3 , y 4 ) space is particularly interesting. In FigJH] and Figj9] we 
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show distributions for the pion-pion rescattering and double-diffractive contributions, re- 
spectively. In this calculation the cut-off parameter A 2 ,^ E = 1 GeV 2 . The cross section 
for the pion-pion rescattering drops quickly with the center-of-mass energy. The rescattered 
pions are emitted preferentially in different hemispheres, ir + at positive 2/3 and n~ at neg- 
ative 2/4 or 7r + at negative y 3 and 7r~ at positive 2/4- The bare (without absorption effects) 
cross section for the double-diffractive contribution grows with energy. At high energies the 
pions are emitted preferentially in the same hemispheres, i.e. 1/3,1/4 > or 2/3,2/4 < 0. While 
at low energies (PANDA) both contributions (exclusive double diffraction and pion-pion 
rescattering) overlap, at high energies (RHIC, Tevatron, LHC) they are well separated, i.e. 
can, at least in principle, be measured. 




FIG. 8: Differential cross section in (2/3,2/4) for the pion-pion rescattering contribution for different 
incident energies: W = 5.5 (PANDA), 200 (RHIC), I960 (Tevatron), 14000 (LHC) GeV. 




FIG. 9: Differential cross section in (2/3,2/4) for the double-diffractive contribution for different 
incident energies: W = 5.5 (PANDA), 200 (RHIC), 1960 (Tevatron), 14000 (LHC) GeV. 

At high energies the diffractive contribution seems more interesting. The camel-like shape 
of the j/3,2/4 distribution requires a separate discussion. In our calculation we include both 
pomeron and reggeon exchanges. In Fig{TD] we show the cross section in y n = 2/3 = 2/4 for all 
ingredients included (thick solid line) and when only pomeron exchanges are included (long 
dashed line), separately for pomeron- reggeon and reggeon-pomeron exchanges (dotted lines) 
and when only reggeon exchanges are included (dashed line). In this calculation the cut-off 
parameter A^ £ = 1 GeV 2 . At low energies all individual cross sections when isolated 
are comparable. They strongly interfere leading to increase of the cross section. At higher 
energies each of the "isolated" cross section peak in different region of 2/3 or 2/4. The IP ® IP 
cross section peaks at midrapidities of pions, while IP ® IR and M ® IP at backward and 
forward pion rapidities, respectively. When interfering the three components in the ampli- 
tude produce significant (camel- like) enhancements of the cross section at forward/backward 
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rapidities. It would be desirable to identify the camel- like structure experimentally 6 . At 
even more forward/backward rapidities one may expect single-diffraction contributions (e.g. 
diffractive production of nucleon resonances and their decays) not included in the present 



analysis. This will be discussed elsewhere [28 




FIG. 10: Rapidity distribution of pions (tt + or tt~) for different center-of-mass energies. The 
different lines corresponds to the situation when only some components in the amplitude are 
included. The details are explained in the main text. 



In Fig{TT] we compare distributions of pion rapidities y n for exclusive double diffraction 
and high-M^ pion-pion rescattering at the PANDA, RHIC, Tevatron and LHC energies. 
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FIG. 11: Differential cross section da/dy^ for the double diffractive and high-M^ pion-pion rescat- 
tering contributions at the PANDA, RHIC, Tevatron and LHC energies. The solid lines was ob- 
tained with Alff E = 1 GeV 2 and dashed lines with A 2 ^ E = 0.5 GeV 2 . 

In Fig{T2]we show the two-pion invariant-mass distribution at the PANDA, RHIC, Teva- 
tron and LHC energies. At the lowest energy the pion-pion rescattering and the double- 
diffractive components strongly overlap. While the double-diffractive component dominates 
at low two-pion invariant masses, the pion-pion rescattering components dominates at large 
invariant masses. This dependence can be used to improve purity of one of the two com- 
ponents by imposing extra cuts. At the Tevatron and LHC energies the double-diffractive 
component dominates over the pion-pion rescattering in the whole range of M n7T . 

Finally we show distributions in transverse momentum of the pions. The double- 
diffractive component dominates over the pion-pion rescattering component. The pions 
are produced preferentially back-to-back. The smearing in p t7T around zero as well as with 



6 The ALICE experiment seems to be able to study the dependence because of the much lower threshold 
on pion transverse momenta. 
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FIG. 12: Differential cross section da/dM nn for diffractive and high-M^ pion-pion rescattering 
contributions at the PANDA, RHIC, Tevatron and LHC energies. The solid lines was obtained 
with ^-offE = 1 GeV 2 and dashed lines with A 2 ^ B = 0.5 GeV 2 . 



respect to <p = it (relative azimuthal angle between charged pions) is caused by the disbalance 
of transverse momenta of exchanged pomerons and/or reggeons from both proton/antiproton 
lines. Whether the distributions can be measured at the LHC requires a Monte Carlo studies 
of the ALICE detector. 
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FIG. 13: Differential cross section da/dpt^ for diffractive and pion-pion rescattering contributions 
at the PANDA, RHIC, Tevatron and LHC energies. The solid lines was obtained with A 2 ^^, = 1 
GeV 2 and the dashed lines with A 2 , f F = 0.5 GeV 2 . 



VII. OUTLOOK 



A. Beyond the Born approximation 

In the present, intentionally simplified, analysis we have performed calculation in the 
Born approximation with the form factor parameter roughly adjusted to existing "low- 
energy" experimental data. In a more microscopic approach one has to include higher-order 
diagrams shown in FigfLJ] and in Figfl5l 

The first type of the interaction was studied e.g. for three-body reactions. For the 
four-body reaction discussed here a similar effect is expected, i.e. large energy-dependent 
damping of the cross section which is often embodied in the soft gap survival probability. 

When going from the Born (Figj3j) to the diagrams with the pion-pion FSI (Fig JT5|) the 
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FIG. 14: Diagrams representing the absorption effects due to proton-proton interaction. 




FIG. 15: Diagrams representing pion-pion final state interaction. 



following replacement is formally required: 
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ff— shell 



(k 3 k 4 



PsPi) 



(7.1) 

where the sum runs over different isospin combinations of pions. In general the integral above 
is complicated (singularities, unknown elements), the vertex form factors (A and B) with two 
pions being off-mass-shell are not well known, and even the off-shell matrix element is not 
fully under control. Usually a serious simplifications are done to make the calculation useful 
on a practical level. Limiting to the S-w&ve (L = 0) one can correct the Born amplitude by 
a phenomenological function which causes an enhancement close to the two-pion threshold 
and damping at M nn ~ 0.8 GeV. Dealing with higher partial waves is more complicated. At 
even larger M nn the interaction becomes absorptive and was not much studied. Some work 



can be found in Ref. [lOj. Clearly much more theoretical afford is required. 

The second type of diagrams leads approximately to a redistribution of the strenght but 
seems to modify the pion-pion integrated cross section very little 0]. The effect of pion- 
pion FSI must be, however, included if the spectrum of invariant mass is studied. At high 
invariant masses one may expect also a strong damping due to absorption in the pion-pion 
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subsystem. Only low- invariant-mass spectra were studied in the past experiments 27j • The 
experiments at LHC could study the potential damping of large-mass dipion production and 
therefore could shed more light on the not fully understood problem of absorption effects in 
a few-body hadronic systems, so important in understanding e.g. the exclusive production 
of the Higgs boson discussed recently in the literature. 



B. Other not included processes at high energies 

Up to now we have discussed only central double- diffract ive (CDD) contribution to the 
pp —* pp7T + 7T~ reaction. In general, there are also contributions with diffractive single or 
double proton/antiproton excitations followed by the resonance decays shown in Figs JTB1 and 
[T71 The first mechanism contribute both to the pp — > pptt + tt~ and pp — > ppn + ir~ reaction 
while the second mechanism only to the pp — > ppn + TT~ reaction at high energy 7 . 
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FIG. 16: Resonance contributions leading to the pp — > ppir + TT channel through diffractive single 
resonance excitation (DSRE). 
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FIG. 17: Resonance contributions leading to the pp — ► ppix^ix channel through diffractive double 
resonance excitation (DDRE). 

Can these processes be separated from the CDD contribution. The general situation at 
high energy is sketched in FigJT8l The discussed in this paper CDD contributions lays along 
the diagonal = y± and the classical DPE in the center y 3 w y 4 . The diffractive single 
resonance excitation (DSRE) contribution 8 is expected at ^3,3/4 ~ ybeam or 2/3,2/4 ~ ytarget, 
i.e. situated at the end points of the CDD contribution. The diffractive double resonance 
excitation (DDRE) contribution is expected at (y 3 ~ y beam and y 4 ~ ytarget) or (y 3 ~ ytarget 
and 1/4 ~ ybeam), i-e. well separated from the CDD contribution discussed in the present 
paper. The Tevatron is the only place where one could look at the DDRE contribution, 
never studied so far at high energies, when it is clearly separated from other mechanisms 
(CDD, DSRE). 



7 At low energy the double A isobar excitation contribute to pp 

8 The Roper resonance excitation is a good example. 
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FIG. 18: A schematic localization of different mechanisms for the pp — > ppir + 7r or pp —* ppir + ir 
reactions at high energies. The acronyms used in the figure are explained in the main text. 

VIII. CONCLUSIONS 

We have calculated several differential observables for the exclusive pp —>■ ppir + 7r~ and 
pp —>■ pp7T + 7T~ reactions. Both double diffractive and pion-pion rescattering proceses were 
considered. The full amplitude was parametrized in terms of subsystem amplitudes. Only 
continuum processes (classical DPE) were included in the present analysis. 

In the first case the energy dependence of the amplitudes of irN subsystems was 
parametrized in the Regge form which describes total and elastic cross section for ttN scat- 
tering. This parametrization includes both leading pomeron trajectory as well as subleading 
reggeon exchanges. Even at relatively high energies the inclusion of reggeon exchanges is 
crucial as amplitudes with different combination of exchanges interfere or/and irN subsys- 
tem energies can be relatively small W n N < 10 GeV. The latter happens when y % +,y w - 3> 
or y^+iU^- <C 0. In this region of the phase space one can expect a competition of sin- 
gle diffractive mechanism. In the literature mainly total single diffractive was calculated. 
We leave the estimation of the SD mechanism contributions to the pp-K + n~ channel for a 
separate study. 

The integrated cross section of the central double- diffractive component grows slowly 
with incident energy if absorption effects are ignored. In principle, the absorption effects 
may even reverse the trend. 

In the second case the pion-pion amplitude was parametrized using a recent phase shift 
analysis at the low pion-pion energies and a Regge form of the continuum obtained by 
assumption of Regge factorization. The factorization assumption is made to estimate the 
process contribution. 

The two contributions occupy slightly different parts of the phase space, have different 
energy dependence and in principle can be resolved experimentally. The interference of 
amplitudes of the both processes is almost negligible. 

The energy dependence of the "diffractive" central production of two-pions is quite dif- 
ferent than the one for elastic scattering, single- or double-diffraction. This is due to the 
specificity of the reaction, where rather the subsystem energies dictate the energy depen- 
dence of the process. 
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At high energies we find a preference for the same hemisphere (same-sign rapidity) emis- 
sion of 7r + and 7r~. At ISR energies the same size emission is about 50 % while at LHC 
energies the same hemisphere emission constitutes about 90 % of all cases. 

In the present analysis we have excluded several resonance contributions. Formally they 
belong to the category (c) and not (d) which we consider in the present analysis. But the 
distinction between the different categories is a bit arbitrary and may be quite involved 
experimentally. Further work is required to estimate contribution of such a process. This 
clearly goes beyond the scope of the present analysis but will be done in the future. 
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